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ABSTRACT 

The asymptotic response of donor stars in interacting binary systems to very 
rapid mass loss is characterized by adiabatic expansion throughout their interiors. 
In this limit, energy generation and heat flow through the stellar interior can be 
neglected. We model this response by constructing model sequences, beginning 
Ph' with a donor star filling its Roche lobe at an arbitrary point in its evolution, 

O ■ holding its specific entropy and composition profiles fixed as mass is removed 

'^ ' from the surface. The stellar interior remains in hydrostatic equilibrium. Lumi- 

^ ■ nosity profiles in these adiabatic models of mass-losing stars can be reconstructed 

from the specific entropy profiles and their gradients. These approximations are 
>- ' validated by comparison with time-dependent binary mass transfer calculations. 

We describe how adiabatic mass loss sequences can be used to quantify threshold 






O ■ conditions for dynamical time scale mass transfer, and to establish the range of 

post-common envelope binaries that are allowed energetically. 
Q ■ In dynamical time scale mass transfer, the adiabatic response of the donor 

Q ■ star drives it to expand beyond its Roche lobe, leading to runaway mass transfer 

and the formation of a common envelope with its companion star. For donor stars 
^ ■ with surface convection zones of any significant depth, this runaway condition 

'^ ■ is encountered early in mass transfer, if at all; but for main sequence stars with 
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radiative envelopes, it may be encountered after a prolonged phase of thermal 
time scale mass transfer, a so-called delayed dynamical instability. We identify 
the critical binary mass ratio for the onset of dynamical time scale mass transfer 
as that ratio for which the adiabatic response of the donor star radius to mass 
loss matches that of its Roche lobe at some point during mass transfer; if the 
ratio of donor to accretor masses exceeds this critical value, dynamical time scale 
mass transfer ensues. 

In common envelope evolution, the dissipation of orbital energy of the binary 
provides the energy to eject the common envelope; the energy budget for this 
process consists essentially of the initial orbital energy of the binary and the 
initial self-energies of the binary components. We emphasize that, because stellar 
core and envelope contribute mutually to each other's gravitational potential 
energy, proper evaluation of the total energy of a star requires integration over 
the entire stellar interior, not the ejected envelope alone as commonly assumed. 
We show that the change in total energy of the donor star, as a function of its 
remaining mass along an adiabatic mass-loss sequence, can be calculated either 
by integration over initial and final models, or by a path integral along the 
mass-loss sequence. That change in total energy of the donor star, combined 
with the requirement that both remnant donor and its companion star fit within 
their respective Roche lobes, then circumscribes energetically possible survivors 
of common envelope evolution. 

Subject headings: binaries: close — stars: evolution — stars: interiors — stars: 
mass loss 



1. Introduction 

Mass transfer is the signature feature of close binary evolution. In relatively wide 
binaries, this mass transfer can occur via capture of some part of the stellar wind of 
the donor by the accretor, but in close systems, the dominant mode of mass transfer is 
by Roche lobe overflow (RLOF). Hydrodynam ic estimates of Roche lobe overflow rates 



(jPaczyhski Sz Sienkiewica 1 19721 : lEggletonl |2006| ) show that mass transfer rates are typically 



of order 



--l!(t) • 



where Mi is the mass of the donor star, Porb the orbital period of the binary, i?L the effective 
Roche lobe radius, which is the radius of a sphere of the same volume as the interior of the 
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Roche lobe, ARi = Ri — Ri^ the donor's radius -Ri minus its Roche lobe radius -Rl, and 
n the donor's effective polytropic index. This mass transfer rate is exquisitely sensitive to 
ARi] so long as mass transfer occurs on a time scale much longer than the dynamical time 
scale of the binary, Ri ^ Ri^. 

We can examine the stability of mass transfer by comparing the responses of donor star 
radius and Roche lobe radius to perturbations in the mass transfer rate. So long as ARi ^ 
i?L, hydrostatic equilibrium prevails through the interior of the donor star, and dynamical 
terms become important only in the vicinity of the inner Lagrangian point. Variations in 
the radius of the donor can be separated into mass- and time-dependent parts: 

The first term on the right-hand side represents not only the evolutionary expansion (or 

contraction) of the donor star in the absence of mass transfer, but also the expansion or 

contraction due to thermal relaxation if, for example, mass transfer has driven the donor out 

of thermal equilibrium. The second term on the right-hand side represents the hydrostatic 

response of the radius to changing mass, with thermal or chemical diffusion and nuclear 

evolution suppressed; that is, it represents the adiabatic response of the donor star to mass 

loss. We will henceforth label this logarithmic derivative by the subscript 'ad' instead of t, 

and define 

dlnRi 



The Roche lobe radius of the donor depends on three quantities (assuming a circular orbit) 
- its mass Mi, the total mass of the binary M, and the orbital angular momentum of the 
binary J: 

where ri^{Mi/M) is the Roche lobe radius in units of the binary separation, and is a function 
only of the binary mass ratio, or as written here, of the ratio of donor to total mass. We can 
thus write 

^i^^^-fH^) ^i-^^+(ll^) ^i-^+(^) ^1^^- (^) 

In conservative mass transfer, as that term is commonly used, M and J are conserved 
quantities, and therefore the latter two terms on the right-hand side of Eq. ([5]) vanish. 
However, more generally, a binary may suffer background mass loss through a stellar wind 
from the companion star, or angular momentum losses through, for example, stellar winds or 
gravitational wave radiation. Moreover, systemic mass and angular momentum losses may 
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contain contributions arising from a loss of part of the mass transfer stream between binary 
components; we can call these consequential losses. We therefore decompose the terms in 
d In M and d In J into terms linearly independent of d In Mi and given the subscript Mi , and 
those consequential losses linearly proportional to d In Mi and given subscript t, by analogy 
with Eq. (|2]). We then write 



dlni?L= — TT^] dt+ — -— ^ dlnMi 



dt 



where 



(91ni?L 
dt 



Ml 



d\nM 



Ml, J 



Ml 



d\nM 



a In Ml 



(6) 



dt 



+ 



Ml 



(91ni?L 
d\nJ 



and 



aini?L 
a In Ml/, 



91n/?L 
aini?L 



91n/?L 

91nM 

(9 In J 



Ml, J 



•9 In J y^,. M V-^lnMi^' 



By analogy with Eq. ([3]) above, we define 



Cl = 



a In Ml 



Afi,A/ 



91nM 
9 In Ml 



91nJ 
9t 



(7) 



Ml 



(9) 



Combining Eqs. 02]) and ([6]), we obtain 



dlni?i/i?L 



dlnRi 



Ml 



91ni?L 
9t 



Ml 



dt + (Cad-CL)dlnMi 



(10) 



Now suppose that the donor star exactly fills its Roche lobe. In order for mass transfer 
to occur, Ri must increase faster (or decrease slower) than i?L would if it were not for mass 
transfer; that is, the coefficient of dt in Eq. (fTOl) must be positive. If the mass transfer rate is 
self-limiting, the stellar radius will parallel the Roche lobe radius as mass transfer proceeds, 
i.e., dln(i?i/i?L) ~ 0. This condition then requires that the second term on the right-hand 
side of Eq. fITOl) be negative. Since din Mi < 0, it follows that self-limiting (stable) mass 
transfer can only occur if (ad > Cl- In that case, we can solve directly for the mass transfer 
rate: 



Ml 



Ml 
Cad — Cl 



dluRi 



Ml 



(91ni?L 
dt 



Ml 



:ii) 
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Typically, Cad and (j^ are of order unity, and mass transfer is driven on a time scale Ri/Ri. 
If the donor remains in thermal equilibrium, mass transfer proceeds on a nuclear time scale. 
If, on the other hand, the donor has been driven out of thermal equilibrium by mass loss, 
it typically expands much more rapidly, on its thermal time scale, and accordingly mass 
transfer proceeds on that thermal time scale. In some circumstances, such as among the 
shorter-period cataclysmic variables or double-degenerate binaries, the driving time scale 
may be that characterizing contraction of the Roche lobe by angular momentum losses. 

Dynamical time scale mass transfer occurs if Cad ^ Cl- Both right-hand terms in Eq. fllOl) 
are then positive, and the donor star radius is then driven by simple hydrostatic relaxation 
far beyond its Roche lobe radius. In this case, the mass transfer rate can in principle 
accelerat e up to a truly dynamical rate (cf. Eg. ([I| ) . Thr ee-dimensional simulations of this 
process fJRasio &: Livid Il996l : iRicker &: Taanj |2006| . 120081 ) show that the donor envelope is 
tidally disrupted, engulfing the companion star while tidal torques from the lagging, distorted 
envelope plunge that companion and the core of the donor i nto orbits deepl y embedded in 
that envelope. The binary enters common envelope evolution (jPaczyhskilll976l ) . The common 
envelope itself lags in rotation behind the embedded cores, which spiral inward as they lose 
orbital energy and angular momentum to the envelope. Ultimately, either the dissipation of 
orbital energy suffices to eject the common envelope, revealing the core of the original donor 
and its companion staio, or the cores merge within the common envelope. Examples of the 
products of common envelope evolution include planetary nebulae with close binary nuclei, 
cataclysmic variables, and other short-period binaries containing at least one degenerate 
component. It is important to recognize that common envelope evolution must be very 
rapid, lest radiative losses quench common envelope ejection, and that energy sources other 
than orbital e nergy dissipatio n are unlikely to contribute significantly to the energy budget 
of the system (IWebbinkll2007l ). Adiabatic mass loss sequences can therefore define the energy 
difference between initial and final states. This energy difference must be supplied by orbital 
energy dissipatioiu. 

Useful insights into the be havior of donor stars under ad iabatic mass loss can be gleaned 
from very simplified models. iHjellming fc WebbinkI (119871 ) investigated the properties of 
polytropic models with power-law equations of state, including complete polytropes, com- 



^ Because the common envelope has a much higher specific entropy than the surface of the accreting 
companion, it is much hotter and more tenuou s at that interface. The c ompanion can accrete very httle 
matter during common envelope evolution 



see 



Hiellming fc TaamI (|l99lh . 



^common envelope sys t ems must, of course, also observe angular momentum conservation (cf. 
Nelemans et al.l 12000 . I2OOII : iNelemans fc Tout! 120051 ) . but this is almost always a less stringent constraint 



on the final system. 
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3 cores with n 



I envelopes and 7 = | 
envelopes with point mass cores). 



equations of state), and 
These models are rele- 



posite polytropes (n 

condensed polytropes (n = | 

vant for some main-sequence stars and giant-branch stars, and illuminate the qualitative 

(and, within lira i ts, qu antitative) behavior of those stars. They were further employed by 

Soberman et al.l ( 119971 ) to explore stability criteria for binary mass transfer. However, these 

models leave many important evolutionary stages of a donor unexplore d, and they leave 



considerable room for improvement even where they may be applic able (JEggleton fc Tout 



19891 : iHan et all 120021 : IPodsiadlowski et al.ll2002l : IChen. k Hani 120081 ). A more so phisticated 
approa ch, using realistic stellar models, is both desirable and possible, as shown by lHjellming 
( ll989al Jbl). This paper represents a step toward filling this need. 



We describe here a method of construction of stellar adiabatic mass loss sequences. In 
a given sequence, the initial model correspo nds to a fu l l stella r model at an arbitrary state 
in its evolution. As in the pioneering work of lHjellmingI (jl989aU bl). mass is removed from the 
surface while stellar interior at each point relaxes adiabatically to a new state of hydrostatic 
equilibrium. In Section 2, we describe the modifications to the normal set of stellar structure 
equations that are needed to construct adiabatic mass loss sequences; how one can recover the 
internal luminosity profile, and hence thermal relaxation rate, from the adiabatic sequences; 
and how we have implemented these concepts in a numerical code. Adiabatic mass loss 
sequences constructed in this way yield Cad, fixing the critical mass ratio for dynamical time 
scale mass transfer, but they also reveal thresholds for a delayed dynamical instability arising 
when the donor is stripped down to a nearly isentropic core. These applications are described 
in Section 3. In Section 4, we discuss the calculation of the total energy of a donor star as it 
loses mass adiabatically, and the application of those results to determining in what orbit, 
if any, a binary can survive common envelope evolution. We conclude with a brief appraisal 
of the relevance of our results to population synthesis studies in Section 5. 



2. Adiabatic Mass Loss 



Adiabatic mass loss models rest on two basic approximations: that in the limit of rapid 
mass loss, the response of the donor star is everywhere locally adiabatic; and that hydrostatic 
equilibrium nevertheless prevails throughout the bulk of the interior of the donor star. The 
first of these approximations implies not only that heat flow can be ignored, but that (under 
the Schwarzschild criterion for convection) convective boundaries are preserved with respect 
to mass. Both the specific entropy profile and the composition profile of the donor interior 
are therefore fixed along an adiabatic mass loss sequence. The second approximation implies 
that fluid flow velocities are negligible (highly subsonic) through the interior of the donor. 
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Solutions for the outflow of gas from a lobe-filling star show th at this assumpti on breaks 
down only in the vicinity of the inner Lagrangian point (see, e.g., lEggletonl ( 20061 )). 



To illustrate the strength and limitations of the adiabatic approximation, we show in 
Fig. [T]the specific entropy profiles of the donor star at five epochs during mass transfer for a 
binary with initial mass Mu = 2Mq, initial mass ratio (donor/accretor) qi = 3.1, and initial 
orbital period Pi = 0.4 d. As mass loss progresses from epoch a to epoch e, the mass transfer 
rate rises from about 10~^ Mq yr~^ (a thermal timescale) to about 10~^ Mq yr~^ (a quasi- 
dynamical timescale). We see that the entropy profile of the donor remains nearly unchanged 
during mass transfer on quasi-dynamical time scales (from curve c to curve e), but that even 
for the thermal timescale mass transfer, thermal relaxation is largely concentrated in the 
outermost layers, which are immediately stripped away. Little heat transport occurs in the 
deep interior. This pattern of thermal relaxation typifies stars with radiative envelopes, with 
rising entropy driven by energy absorption in the outer envelope (as opacity rises rapidly 
with decreasing temperature during decompression), accompanied by falling entropy in the 
deep interior (as thermal energy generation supplants nuclear sources as they are frozen out 
by decompression). Among stars with deep surface convection zones, the entire convective 
envelope slowly declines in specific entropy, as radiative energy losses from the surface exceed 
the input of energy at the base of the convection zone (these being choked off by the local rise 
in opacity as the underlying radiative zone undergoes decompression). Thus, the adiabatic 
approximation not only defines the asymptotic response of a donor star in the limit of faster- 
than-thermal time scale mass transfer, but serves as a good approximation to the deep 
interior response even during thermal time scale mass transfer. 



2.1. The Equations of Stellar Structure 

In principle, we can draw the initial model for an adiabatic mass loss sequence from 

any point in the evolution of the donor star. Following common practice in close binary 

evolutionary models, we neglect rotational and tidal effects on the structure of that star, 

and treat it as spherically symmetric. The initial model then satisfies the usual set of four 

first-order ordinary differential equations of stellar structure, where four structure variables, 

i.e., pressure P, radius r, temperature T and luminosity L, vary with an independent variable, 

mass m, 

d In P Gm 



dm 47rr^P 

dlnr 1 

dm Anr^p 



(12) 
(13) 



dlnT _ dlnP 
dm dm ' 

., „ f Vr = . „ '^^KJ^ ^A Radiative zone ,,^, 

with V = < levracGmT^ , (15) 

(^ Vconv Convective zone 

Here density p, opacity k, adiabatic temperature gradient Va = (91nT/91nP)^, nuclear 
energy generation rate e, neutrino loss rate e^, and specific heat at constant pressure Cp = 
{ds/d\nT)p are functions of P, T, and the abundances of various nuclear species X,. Sub- 
script s means entropy. G, c and a are of course the universal gravitational constant, the 
speed of light and the radiation constant, respectively. Vr is temperature gradient that 
would prevail in the absence of convection, while the temp erature gradieii t in a convective 



zone can be written in a common mixing- length treatment ( lEggletonlll983l ) as 



V = V ^^^°°^ (U) 

V conv ^f I ^ ■ y^l J 

The characteristic convective velocity f conv is the real root of the cubic equation 

^^^conv + ^conv + ^^conv = a^y.CpT{y, - V .) , (18) 

X OiHp 

where Hp is the pressure scale height, a the dimensionless mixing- length parameter, and 

AacT^ , , 

Four boundary conditions are required to close this set of equations. At the center (m = 0), 

r = , (20) 

L = ; (21) 

and at the surface {m = M), 

L = Anr^aT^ , (22) 

^ (-'gas + 2-* radj _ 2 , , 

9 3 ' 

where the latter surface boundary condition approximates a classical gray atmosphere. The 
four equations are solved for a given distribution of Xj(m), which is updated according to 
nuclear reactions and convective mixing. Stellar mass loss can be included by applying the 
surface boundary conditions at a moving surface such that dM = M dt. Roche lobe overfiow 
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in a binary system can likewise be accommodated either by expressing M explicitly in terms 
of the degree of Roche lobe overflow (see, for example, the Appendix to this paper), or (in 
the approximation that R = Ri,) by rewriting the structure equations above in terms of r 
rather than m as the independent variable, and applying the surface boundary conditions at 
r = Ri^im). 

In adiabatic mass loss models, the equations of energy transport and conservation 
(Eqs. (ini) and (TTB]) . respectively) are supplanted by algebraic conditions fixing the specific 
entropy and composition profiles: 

s{m) = So{m) (24) 

X,(m)=X,,o(m) , (25) 

where the zero subscripts denote the values of these quantities at the same mass coordinate 
in the initial model of the sequence. Eqs. (fT2|) . (fT3l) . (l24|) and (1251) . together with appropriate 
boundary conditions (identified below), are sufficient to describe completely the state of the 
gas at any point in the interior of the donor. The local temperature, T = T{P, s, Xj), follows 
implicitly. 

To understand the internal relaxation processes set in motion by adiabatic mass loss, 
we would also like to recover the luminosity profile through the interior of the donor star. 
This we can accomplish by writing the local energy fiux through the stellar interior in terms 
of the local entropy gradient, which is a known quantity. Given the local state of the gas, 
the entropy gradient can be written in terms of the pressure and composition gradients: 



ds 
dm 



ds 
dP 



T,Xi 



dP 

dm 



+ 



ds 



d\nP 



Cr 



V 



CpfV- V, 



+ 

T,X, 

ds 
d\nP 

dlnP 

dm 



ds\ dT sr^ f ds 
'Wr I _ _ H^ 2-^ \dX. 

dlnP 



dT ^ 
ds 



P,X, 



d\nT 



T,X, 



P,X, 

ds 



dlnT 
dlnP 



+Ei 



d\nT 

/ ds 
[dX.. 



-1 



P,X^ 



P,T 



dm. 

dlnP 

dm, 

dX 



dX, 
prp dm, 



P,T 



El 



/ ds 
[dX, 



P,T 



dXj 

dm 

dX, 
dm, 



dm. 



(26) 



where we have invoked the definitions of specific heat, Cp = {ds/d\nT)p^Xi, and of the adi- 
abatic temperature gradient, Va = {d\nT/d\iaP)s^Xi = —{ds/dhiP)T,xSds/d\iiiT)p^Xi)~^- 
As overlying layers are stripped away, successively deeper mass layers are brought toward 
the stellar surface, pressure in those layers decreases, and |dlnP/dr7i| increases accordingly. 
The ambient temperature gradient, V, is therefore driven toward the adiabatic temperature 
gradient, Va, since ds/dm, and dXj/dm are fixed in Equation fl26l) . 
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Let us define an effective entropy gradient, ds/dm, as 

dm dm ^—^ \dXi) pj, dm ^ dm 

Then the ambient dimensionless temperature gradient, V, can be written as 

^ ^ dS 1 /dlnPV^ ^ 4:TTr^P ds 

V = VaH = Va . (28) 

dm Cp \ dm J GmCp dm 

The value of ds/dm, can readily be obtained at each layer of an adiabatic mass loss model, 
and so V is a known quantity. If ds/dm > 0, then V < Va; such a layer is stable against 
convection. In this case, Vr = V, and 

IGiracGmT^^ , . 

If ds/dm < 0, on the other hand, then V > Va. Such a layer is convectively unstable 
and energy is transported by both radiation and convection. The temperature gradient that 
would prevail in the absence of convection is (cf. Eqs. [T71 [181 and lTOj) : 

^ 3 

27H],a'^CpTV 
where 



'^+ ..rr.L rr.^ ^\ (30) 



and 

^ ( dg/dm \ (TIHWCpTVA 

\Cp dlnP/dmy' V 4^2 ; ■ ^ ' 

The total luminosity is then 

_ levracGmT^ 
L - ^^ Vr , [ii] 

where Vr is calculated from Eqs. ( 130|) . (!3T!) and (!32l) . Thus the additional information needed 
to reconstruct the luminosity profile in an adiabatic model is already carried in the entropy 
and composition gradients. 

It remains to specify the boundary conditions for the adiabatic mass loss models. At 
the center of the star (m = 0), clearly 

r = (34) 

as before. At the surface of the star, however, there is a fundamental inconsistency between 
the assumption of adiabatic relaxation (implying no radiative losses) and the finite effective 
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temperature that finite photospheric optical depth inevitably demands. In the outer atmo- 
sphere of a star, at small optical depth, radiative relaxation always proceeds more rapidly 
than dynamical relaxation, so that purely adiabatic relaxation is impossible. However, at 
large optical depth, the diffusion approximation remains valid, and it fixes the stellar lumi- 
nosity in terms of the ambient entropy gradient, as we have seen above. We therefore resolve 
this inconsistency by insisting on continuity of luminosity at the stellar surface. Substituting 
the blackbody law into the radiative diffusion equation (Eq. (133|) ). we obtain for the surface 
boundary condition 

(35) 



_ 3 Pk 

Vr = . 

16 g 



This condition is tantamount to one relating the entropy gradient to the entropy itself at 
the photosphere. 



2.2. Numerical Implementation 



The numerical code to solve the equations of the adiabatic mass loss mo del in sect i on [27 



is written in FO RTRAN95, ba. sed on the stellar evolution code developed by lEggletonl (11971 



19721 . 119731 ) and iPaxtonI (120041 ). This code has b een widely used and is upd a ted w i th th e 
latest input physics, as described by iHan et al.l (11994 120031 ) and iPols et al.l (119951 Il998l ). 
The essential features of this code have been retained in our adiabatic mass loss model, 
namely (a) the use of an adaptive, moving, non-Lagrangian mesh; (b) the treatment of 
both convective and semiconvective mixing as diffusion processes; and (c) the simultaneous, 
implicit solution of both the stellar structure equations and (in the evolutionary models) the 
chemical composition equations, including convective mixing. 



Our code uses an equation of state, described by lEggleton et al.l (Il973[ ) and iPols et al. 



( 119951 ). that includes pressure ionization. Coulomb interactions and disso c iation of molec- 



ular hydrogen. The opacity tables are taken from lAlexander fc FergusonI (Il994a 



in the 
in the 



low-temperature limit (where molecules become important), from lltoh et al.l ( 1198 
high- density, high-temp e ratur e regime (where electron conduction becomes important), and 
from JRogers fc Iglesiad (119921) elsewhere. Th e nucl ear r eaction rates ar e inter polated in 
temperature, mostly from ICaughlan fc Fowlerl ( 1l988l ) and ICaughlan et al.l (Il985l). Electron 



screening theory, including strong and intermediate screening, comes from iGraboske et al. 



( 119731 ). The neutrino energy loss rates resulting from pair, photo-, pla sma, bremsstrahlung, 
and recombination neutrino processes originate from lltoh et al.l (Il996[ ). 



The adaptive, non-Lagrangian mesh automatically dist ributes rnesh p oints at uniform 
intervals in a function, Q, of the local structure variables (lEggletonlll97ll ). The choice of 
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function Q is in principle arbitrary, subject only to the constraint that it vary monoton- 
ically from center to surface of the stellar model. Its introduction involves rewriting the 
stellar structure equations in terms of fc, the mesh point number, as the independent (radial) 
variable, and introducing an additional 'structure' equation for the mesh function: 

dg Q(iv)-Q(i) ^ ^ .„„. 

—— = = constant , (do) 

ak N — 1 

where N is the number of mesh points in the model. In our case, we employ 400 mesh points, 
distributed according to the mesh function 



Q = 0.05 ln(P/dyn cm"^) + 0.15 In [(P/dyn cm"^) + 10^^] + 0.45 In [T/{T + 2 ■ lO^K) 

+ In [0.02 m2/7(0.02 m^/^ + m^/^)] - 0.05 In [l + (r/10^ cm)^] , (37) 

where rUc = 3.5p(P/G'p^)'^/^. We find that his mesh function provides good resolution of 
structural features of interest throughout our evolutionary and mass loss sequences. Taking 
the derivative of Eq. f l37|l . and substituting Eqs. f lT2| [T3| andlH]), with V determined from 
Eq. ( 125]) in lieu of Eq. (IT^ . we can find dQ/dm. Using Eq. (15^ . we can then cast the 
structure equations in terms of k rather than m as the independent (radial) variable. 

In practi ce, we use In f ( a deg e neracy pa r amete r related to the electron chemical po- 



tential — see lEggleton et al.l ( 1l973l ): lEggletonl (120061 ) ). and ln(T/K) as our state variables 



and ln(r/10^^ cm), {m/10^^ g), and Q as our global variables. The differential equations 
for P{k), r{k), m{k), and Q{k) are written in difference form, as usual, but we include also 
the algebraic equation so{m{k)) = s{f,T,Xi) in the solution matrix as a convenient tool for 
satisfying the entropy constraint. The two additional differential equations require boundary 
conditions. At the surface {k = 1), 

Vr = ^- , (38) 

16 g 

m = M . (39) 

At the center {k = N), we should have r = 0, as before, and m = 0, but to avoid singular 
behavior in that limit, the central mesh point is offset: 

m = -- , (40) 

dr 

As intimated above, the construction of adiabatic mass loss sequences involves the inter- 
polation in this adaptive mesh of the specific entropy s at the corresponding mass coordinate 
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in the initial model of the sequence. Near the surface of an initial model, s typically varies 
very rapidly with mass, while the mass increment from one mesh point to the next becomes 
a very small fraction of the stellar mass. Because of the limitations of numerical accuracy 
(we typically require that the average correction to structure variables in a model deemed 
'converged' not exceed one part in 10^), the situation arises frequently that the nominal error 
in mass coordinate m is comparable with (or even exceeds) the mass difference between one 
mesh point and the next. In this case, the nominal mass coordinates of mesh points near the 
surface may not be well-ordered, making the accurate interpolation in mass of the specific 
entropy and its gradient impossible. We avoid this problem altogether by using not the mass 
coordinate m, but the mass complement, fn = m — M, where M is the total mass of the star, 
over the outer half of our mesh. Our convergence criterion, now applied to fn, ensures a mesh 
well-ordered in mass. For similar reasons of numerical accuracy, we tabulate ds/dm, derived 
from the local stellar luminosity as described above, for the initial model of each adiabatic 
sequence, and interpolate in this gradient, rather than numerically differentiating s and X 
themselves, in reconstructing the internal luminosity profiles of models in our adiabatic mass 
loss sequences. 

As a demonstration of the ability of our adiabatic mass loss code to reproduce faithfully 
the response of a donor star undergoing rapid mass loss, we have calculated an adiabatic 
mass loss sequence, beginning at model (c) of the time-dependent mass transfer calculation 
shown in Fig. [H At this point, the mass loss rate from the donor has grown to roughly a 
thermal time scale rate. In Fig. [2] we show comparisons of the interior density, luminosity, 
temperature, and radius profiles in the adiabatic models with those in the time-dependent 
calculation when both models have reached point (e) of Fig. [1] Only very minor differences 
appear in layers very close to the surface, but in general these two calculations are in excellent 
agreement, giving strong support to the use of adiabatic mass loss sequences to model rapid 
mass transfer. 



3. Thresholds for Dynamical Instability 

3.1. Physical Significance of Qd 

An important application of adiabatic mass loss sequences is the evaluation of threshold 
conditions for the onset of dynamical time scale mass transfer. As discussed in Section [H this 
condition corresponds formally to Cad < Cl; in that case, a sufficiently rapid perturbation to 
the donor star (i.e., one faster than a thermal time scale) can lead to unbridled growth of the 
donor radius beyond its Roche lobe (cf. Eq. [10]). On the other hand, if Cad > Cl the donor 
will transfer mass on a slower time scale, on a nuclear or angular momentum loss time scale 
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if the donor thermal equihbrium radius-mass relation (Ceq = (i91ni?/(91nM)eq) is sufficiently 
steep (Ceq > Cl); or on a thermal time scale if thermal relaxation tends to drive the donor 
beyond its Roche lobe (Ceq < Cl < Cad)- 

Adiabatic responses of stars to mass loss depend strongly on the variation of specific 
entropy with mass in their outermost layers. This behavior is illustrated in Fig. |3] for lower 
main sequence stars. As one progresses down the main sequence, the surface convection zone 
grows deeper, and, as illustrated here, the initial contraction of stellar radius with decreasing 
mass becomes less and less severe. In the limit of fully convective stars (M; = O.SOMq in 
Fig. [3]), their behavior approximates the R/Ru ~ (M/Mii)~^/^ dependence that follows 
from the homologous expansion of an isentropic gas sphere supported by classical ideal gas 
pressure. In contrast, stars on the upper main sequence, with nearly completely radiative 
envelopes, undergo very rapid contraction with adiabatic mass loss, since their high-entropy 
outer layers occupy a disproportionately large fraction of their stellar volume. Broadly 
speaking, then, stars with convective envelopes tend toward dynamically unstable mass 
transfer, while those with radiative envelopes tend toward stability. 

The situation in reality is somewhat more complex. As is evident in Fig. [31 the radius- 
mass exponents Cad, Ceq, and Cl are not constant over the course of mass transfer, but vary 
continuously, in some cases very rapidly. This is true particularly of Cad in the case of a 
star with a radiative envelope as it first fills its Roche lobe. Typically, the specific entropy 
rises extremely rapidly through the surface layers of such a star (see Fig. [T]), which therefore 
decrease rapidly in density. Loss of even a small amount of mass leads to rapid reduction 
in the stellar radius, so Cad is initially very large and positive. Donor stars with radiative 
envelopes are therefore typically strongly stable against dynamical time scale mass transfer, 
at least initially. However, as mass loss proceeds, and the surface layers are stripped away, 
entropy gradients flatten, and the initially rapid stellar contraction moderates considerably. 
If the deep interior lacks strong entropy gradients (a circumstance typically prevailing during 
core hydrogen burning), then as the donor star is stripped down to this region it may begin 
expanding with further mass loss, much as an (isentropic) convective star expands with 
decreasing mass. Under these circumstances, the donor star may belatedly become unstable 
to delayed dynamical mass transfer. 



3.2. The Delayed Dynamical Instability 

The adiabatic response of a 3.2 Mq terminal main sequence star, shown in Fig. HI 
illustrates the nature of this delayed dynamical instability. This star has a deep radiative 
envelope, and we see that in the adiabatic limit, its initial response is indeed extremely rapid 
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contraction. Recalling that even in thermal time scale mass transfer, little thermal relaxation 
occurs in the deep interior of a donor star (see section [2], Fig. [1]), we can regard the adiabatic 
mass loss sequence as approximating the minimum radius of the donor star as a function 
of its mass. We see that, in this example, a binary with an initial mass ratio gi = 3.0, 
assumed to transfer mass conservatively, never reaches this minimum radius in the course of 
mass transfer; in this case, the mass transfer rate need never exceed a thermal rate in order 
for the donor star to track its Roche lobe through the mass transfer episode. In contrast, 
a binary with an initial mass ratio gj = 4.5 undergoes such rapid Roche lobe contraction 
(again in the conservative approximation) that, despite the rapid initial adiabatic contraction 
of the donor, the Roche lobe radius reaches the minimum radius set by the adiabatic mass 
loss sequence (intersection of the dashed line with the solid line in Fig. H]) when the donor has 
been reduced to a fraction 0.86 of its initial mass. In this case, the system may transfer mass 
at a thermal rate up to that point, but the donor is then no longer able to remain within 
its Roche lobe, and dynamical time scale mass transfer must ensue. This is the delayed 
dynamical instability. At an initial mass ratio gi = 3.757 (dash-dotted curve), the Roche 
lobe radius is tangent to the adiabatic mass-radius relation (at M/Mj = 0.653), a condition 
which then marks the threshold for delayed dynamical instability, in the limit that the donor 
star response is truly adiabatic, and that mass transfer is conservative. 

Since the delayed dynamical instability follows a period of slower, thermal time scale 
mass transfer, in which some measure of thermal relaxation inevitably occurs, the question 
arises, how good is the adiabatic approximation for the threshold of delayed dynamical 
instability? We can glean at least a partial answer to this question by comparing the adiabatic 
results in the example above with a family of time-dependent mass transfer models spanning 
the same range of initial mass ratios, and again assuming conservative mass transfer, but now 
with the full set of stellar structure equations. We again assume conservative mass transfer, 
but now require the donor star to exactly fill its Roche lobe throughout mass transfer. The 
results of this exercise are shown in Fig. El showing mass transfer rates as a function of 
donor mass on the same scale as Fig. H] to facilitate comparison. In this figure, the onset 
of delayed dynamical instability is marked by a rapidly accelerating mass transfer rate as a 
function of decreasing donor star mass. We see that the mass transfer rates in systems with 
initial mass ratios gi = 4.5 and qi = 3.757 both indeed accelerate well beyond the thermal 
time scale rate, while the transfer rate for the system with initial mass ratio gi = 3.0 never 
accelerates beyond a thermal rate. The mass loss rates for the unstable systems g^ = 4.5 
and gi = 3.757 peak at donor star masses M/Mi = 0.905 and M/Mi = 0.814, respectively, 
somewhat earlier than predicted by the adiabatic calculations seen in Fig. m indicating that 
the modest thermal relaxation that occurs in the time-dependent calculation drives the donor 
slightly in the direction of instability. In this case, we may overestimate the critical initial 
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mass ratio for delayed dynamical instability by about 10%, which we may then consider a 
mark of the accuracy and reliability with which we can define this threshold from adiabatic 
mass loss sequences alone. 



3.3. Superadiabatic Expansion 

A different sort of complication arises among luminous giant branch stars, in which 
convection becomes very inefficient in their outer envelopes. The ambient density there falls 
so low that an increasing fraction of the stellar luminosity must be carried by radiation, 
notwithstanding the fact that the local medium is unstable to convection. This condition 
drives the local temperature gradient, V, well above the adiabatic gradient, Va (see sub- 
section [2]T]) . The entropy gradient thus becomes negative over a significant fraction of the 
outer envelope; in effect, the outer envelope forms a cool, dense blanket, suppressing the 
higher-entropy convective envelope contained within it. Removal of this blanketing supera- 
diabatic layer leads to an extremely rapid initial expansion of the star. As seen in Fig. El 
this superadiabatic expansion can be quite dramatic in luminous giant stars, but in fact it 
is a feature of any donor star with a surface convection zone. Among less luminous stars 
(generally, those which have not yet reached the giant branch) , the mass and radial extent of 
the superadiabatic layer is small enough that its effect on threshold conditions for sustained 
dynamical time scale mass transfer are negligible. However, superadiabatic expansion is a 
pervasive feature of luminous giant branch stars. 

The behavior seen in Fig. |6] reveals certain limitations to the utility of spherically- 
symmetric adiabatic mass loss models in analyzing the stability of mass transfer. In reality, 
mass is not removed uniformly over the whole surface of the donor, but from an initially 
small region about the inner Lagrangian point, Li. Moreover, sufficiently close to the stel- 
lar photosphere, radiative relaxation inevitably becomes rapid compared with convective 
turnover time scales, so that conditions for adiabatic outflow cannot be satisfled. Indeed, 
for a giant branch star near Eddington luminosities, the thermal time scale of the entire 
envelope may approach its dynamical time scale. Without full three-dimensional models, 
including both advective and radiative energy transport, we cannot be certain how much of 
this superadiabatic expansion is real, and how much an artifact of our simplifying assump- 
tions. Nevertheless, there exists a real physical basis for this behavior, and the adiabatic 
models show that it is at least energetically possible for superadiabatic expansion to trigger 
mass transfer when the donor photosphere still lies well within its Roche lobe. In the con- 
vective envelope, hydrostatic equilibrium describes an average state of the gas, but in fact 
it consists of rising and descending flows, which may differ considerably from each other in 
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specific energy or entropy, and approach sonic velocity where the mean entropy gradient is 
strongly superadiabatic. That kinetic energy is in principal capable of launching gas through 
the saddle point in gravitational potential at Li even as the mean stellar photosphere lies 
well within the donor star Roche lobe. Rising convective flows may be captured by the 
companion star, with the result that mass transfer rates can fluctuate wildly on convective 
turnover time scalesj. 



3.4. Critical Mass Ratios 



In future papers (Ge et al. , in preparation) , we will present results of a survey of adiabatic 
responses of stars throughout the Hertzsprung-Russell diagram. Inasmuch as one of the 
principal objectives of this survey is to establish the conditions under which a lobe-filling 
donor star of arbitrary evolutionary state is stable or unstable to dynamical time scale mass 
transfer, we describe here the issues which arise in characterizing our results for that purpose. 

We idealize mass transfer between binary components as fully conservative, by which 
we mean specifically that the total mass of the binary 



Mtot = Ml + M2 
and orbital angular momentum of the binary 
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are constant. (We assume, of c ourse , that the orbital eccentricity is throughout negligibly 
small.) Then, using lEggletonl 's ( 119831 ) approximation for the Roche lobe radius. 



rh{q) 



0.49g2/3 



0.6g2/3 + in(l + gl/3 



(44) 



where the mass ratio q is defined, as before, as the ratio of donor star mass to accretor mass, 
we can write the Roche lobe mass-radius relation directly in terms of the binary mass ratio: 



Cl = 
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a In Ml 



J,M 



2 ln(l + gi/3) _ gi/3/(i + ^1/3) (l-q 
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Episodic mass transfer driven by sup eradiabatic expansion was first postulated by Bath (jBathl (|1972 , 



19751 ): see also iPapaloizou fc BathI (|l975l) V 
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Given (^ad from the adiabatic mass loss sequence for a donor star of interest, Eq. fH5|) then 
imphcitly defines a corresponding critical mass ratio, gcrit, satisfying the equation 

Cad = CL(gcrit) (46) 

above which a binary containing that donor star is unstable to dynamical time scale mass 
transfer. 

In reality, mass transfer is never truly conservative, and Eq. (145|) is at best an approxi- 
mation to the dependence of Roche lobe radius on mass ratio. In contrast, to the extent that 
our spherically-symmetric models are appropriate, Cad is intrinsic to the donor star alone, 
independent of the mass or nature of the companion star. In this sense, it is the more fun- 
damental parameter by which we should characterize our results. However, as noted above. 
Cad varies continuously over the course of mass loss, and the decrease in Cad with decreasing 
donor mass is crucial to the existence of a delayed dynamical instability. Conditions for 
the onset of that instability therefore depend on the extent and nature of mass and angular 
momentum losses accompanying mass transfer. For the sake of clarity and simplicity, we 
here and henceforth implicitly assume conservative mass transfer in quantifying thresholds 
for dynamical time scale mass transfer, and denote by gcrit, specifically the threshold initial 
mass ratio for dynamical instability, and define Ccrit = (hiQcTit), the corresponding critical 
donor mass-radius exponent marking this same threshold. 

The evaluation of gcrit is comphcated by the fact that the distinction between a prompt 
dynamical instability and a delayed dynamical instability is to some extent artificial. Even 
if the donor star satisfies the inequality for dynamical instability (Cad < Cl) upon filling its 
Roche lobe, that instability will not be manifested until the mass transfer rate exceeds a 
thermal rate, and becomes too rapid for thermal relaxation to damp its growth. We have 
therefore adopted an approach to evaluating threshold conditions for dynamical instability 
that accounts explicitly, if crudely, for the phase of growth in mass transfer rate up to the 
point at which donor star response can be reasonably approximated as adiabatic. 

As a first approximation of the response of the donor star, we suppose that its structure 
can still be approximated by its adiabatic mass loss sequence, notwithstanding thermal 
relaxation. This is the same assumption employed above in illustrating the threshold for a 
delayed dynamical instability, and is subject to the limitations and shortcomings described 
there. In a prompt dynamical instability, on the other hand, this should be a much better 
approximation, as the growth time scale for the mass transfer rate is generally much shorter 
than a thermal time scale. 

The mass loss rate from the donor is a function of the structure of the outer envelope 
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of the donor star, and of the degree to which it overfills its Roche lobe, 

Ml = Mi{ARi) . (47) 

Details of the calculation of Mi as adopted in this study can be found in the Appendix. At 
the start of a mass loss sequence, the donor is assumed to fill its Roche lobe exactly: 

Rii = Rw (48) 

As mass transfer proceeds, we assume that the Roche lobe follows the radius-mass relation 
for conservative mass transfer appropriate to its initial mass ratio, gi = Mii/M2i. Defining 
mass fraction 

we can write 

where clearly /ii = gi/(l + q{). The donor radius Ri follows the adiabatic mass loss sequence, 
by assumption: 

i?l(Mi)=i?ad(/iMii//ii) , (51) 

where Mi = fiMn/fii follows from conservation of total mass. The desired threshold condi- 
tion for dynamical time scale mass transfer then consists in finding that initial mass fraction 
/ii (or, equivalently, initial mass ratio gi) such that the Roche lobe radius, at its deepest pen- 
etration into the envelope of the donor star, suffices to drive mass transfer at just a thermal 
(Kelvin-Helmholtz) rate. For that rate, we adopt 

GMu , , 

Mkh = -^^ , (52) 

where the subscripts "Ii" refer to the values of the corresponding variables at the beginning 
of the mass loss sequence. 

Let us define i?KH(^i) as the radius within the donor star at which 

M(Akh) = Mkh , (53) 

where 

AKH = i?i(Mi)-i?KH(Mi) . (54) 

The algorithm we have employed to evaluate Eq. (!53|) is detailed in the Appendix to this 

paper. i?KH is then implicitly defined by Eqs. (15T]) - (!54|) as a function of Mii, /x, and /ii. We 

likewise define 

_ /dhii?KH\ 
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The our desired solution simultaneously satisfies the equations 

i?KH(Mii, /i, /ii) = Rhift, /ii) (56) 

and 

CKH(Mii,/i,/ii) = CL(/i) • (57) 

Mil is fixed by the choice of initial model, so these constitute two equations in two unknowns 
(/^ and /ii). The solution value for Hi then corresponds to the minimum initial mass fraction 
above which the binary suffers dynamical instability; that for /i specifies the mass fraction 
at which the instability actually occurs. The corresponding critical initial mass ratio is 
Q'crit = /^i/(l ~ A^i)) and we define 

Ccrit = CL(gcrit) (58) 

(see Eq. HS]) as the initial mass-radius exponent characterizing this critical case. If the binary 
is subject to a prompt instability, /ii and fi differ by a negligible amount, but in a delayed 
dynamical instability the difference provides an approximate measure of the extent of mass 
transfer prior to the onset of dynamical instability. 

In practice, we have found that this approach works quite satisfactorily, and is insen- 
sitive to the vagaries of round-off error that can afflict numerical derivatives like Cad and 
Ckh- However, in the case of donors undergoing pronounced superadiabatic expansion, the 
ambiguities described in Section [3731 lead us to characterize our results somewhat differently. 

We conjecture that, when superadiabatic expansion becomes important, the true thresh- 
old for dynamical instability probably lies between two extremes. If we suppose, on the one 
hand, that a donor only initiates mass transfer when it fills its Roche lobe, as assumed above, 
then the rapid superadiabatic expansion that follows presents a most favorable circumstance 
for dynamical instability. The algorithm detailed above then yields estimates for gcrit.min 
and Ccrit,min, the minimum plausible mass ratio, and corresponding initial mass-radius expo- 
nent, unstable to dynamical time scale mass transfer. If superadiabatic expansion is severe 
(Ckh < ~5/3), no simultaneous solution exists for Eqs. (j56|) and fl57|) in the case of conser- 
vative mass transfer. In this case, gcrit.min -^ 0, i.e., all mass ratios are dynamically unstable, 
and we set Ccrit,min equal to the minimum value of ln(i?KH/-Rii)/ln(Mi/Mii) along the mass 
loss sequence. 

At the opposite extreme, as noted in Section 13.31 it is energetically possible for a donor 
star to initiate mass transfer while its photosphere still lies well within its Roche lobe, given a 
sufficiently strong perturbation. This circumstance we consider the least prone to dynamical 
instability, because sustaining mass transfer at a high rate requires rapid contraction of the 
donor Roche lobe, and therefore a relatively extreme mass ratio. We characterize the initial 
rapid superadiabatic expansion by Apxp, the logarithmic difference between the Roche lobe 
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radius and the stellar radius at the onset of mass transfer. To quantify Agxp, we build a 
pseudo-envelope model for the donor star in question by assigning to its entire convective 
envelope a uniform specific entropy equal to that at the base of the convective envelope 
of the true evolutionary model. By construction, this pseudo-model has no superadiabatic 
surface layer. As seen in Fig. [6] below, the difference between true and pseudo-models is 
evident only very near the surface. Then, identifying the radius of this pseudo-model with 
the Roche lobe radius of the donor, we define 

Ae.p = \og{Ru/Rii) , (59) 

where i?i(Mi) is the radius of the pseudo- model as a function of mass. Agxp is an estimate of 
the fraction by which a donor giant branch star may underfill its Roche lobe and yet initiate 
tidal mass transfer to its companion. In turn, the same algorithm described above for finding 
gcrit and Ccrit, applied now to the pseudo-envelope model, yields estimates for gcrit.max and 
Ccrit,max, the maximum plausible mass ratio, and corresponding initial mass-radius exponent, 
stable against dynamical time scale mass transfer. 

Fig. E] illustrates the application of the algorithms described above to the determination 
of threshold conditions for dynamical time scale mass transfer. In this example, the donor 
is a 1 Mq star near the tip of the giant branch. In response to adiabatic mass loss, the 
radius of such a donor (heavy solid line) undergoes a rapid initial expansion that will trigger 
dynamical time scale mass transfer, provided that the mass transfer rate becomes large 
enough to damp thermal relaxation. In this example, the donor has a very extended envelope 
and short thermal time scale, and so it must overfill its Roche lobe by a significant factor 
in order to drive the mass transfer rate up to Mkh; the depth in the donor at which this 
occurs is shown by the thin solid line. The threshold condition for dynamical time scale 
mass transfer then occurs when the Roche lobe radius just reaches -Rkh- In the example 
at hand, however, expansion is extremely rapid that all mass ratios are unstable: Ccrit,min = 
—27.4 and gcrit,min = 0. In this critical case, M reaches Mkh at log(Mi/Mii) = —0.00027, 
log(-RKH/-Rii) = 0.00742; that is, dynamical instability is virtually instantaneous. 

Also shown in Fig. E] is a pseudo-envelope model sequence for the same 1 Mq giant 
branch star (heavy dashed line). The radius in this model at which M = Mkh is shown 
by the thin dashed line. In the limiting case, it is energetically possible for this star to 
sustain Roche lobe overflow by adiabatic expansion when it initially underfills its Roche 
lobe by Aexp- At the threshold for this case the Roche lobe radius (dotted line) again 
touches the thin dashed line but, because the donor star underfills its Roche lobe, more 
dramatic contraction of the Roche lobe is needed to drive M up to Mkh- The corresponding 
tangent point occurs at log(Mi/Mii) = -0.096, log(^KH/i?ii) = 0.007, where Cl = 0.455 
and q = 0.999, corresponding to gcritmax = 1-661 and Ccritmax = 1-928 at Mi = M^. 
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Our conjecture that the threshold initial mass ratio for dynamical time scale mass 
transfer, gcrit, probably lies within the interval gcrit,min < feit < fcit.max, as these limits 
are estimated above (and likewise for the critical initial mass-radius exponent, Ccrit,min ^ 
Ccrit < Ccrit,max) niust remain a matter of speculation at present. It should be noted that 
the example illustrated in Fig. |6] is a relatively extreme one, and that given the very short 
thermal time scale of the donor, even dynamically-stable mass transfer on a thermal time 
scale could trigger common envelope evolution. 



4. Common Envelope Evolution 

The family of binary stars contains many examples of non-interacting systems con- 
taining a compact component (white dwarfs, neutron stars, or black holes) in an orbit so 
compact that the immediate progenitor of the compact component could not possibly have 
fit within its Roche lobe. These systems must have originated through common envelope 



evolution (jPaczyhskil Il976l ) . Examples include planetary nebulae with close binary nuclei, 
pre-cataclysmic white dwarf-main sequence star pairs, close double white dwarfs, and many 
binary pulsars. These systems are presumably the progenitors to interacting binaries contain- 
ing compact components, either through evolution of a non-degenerate stellar component, or 
through angular momentum loss from the binary via a magnetic stellar wind or gravitational 
radiation, for example. 

Modeling common envelope evolution poses one of the most demanding challenges in 
computational stellar astrophysics. The characteristic dimensions of the common envelope 
may easily exceed those of a compact core embedded within it by factors of 10^ or more, 
corresponding to a ratio of dynamical time scales (envelope versus core) exceeding 10^. Fur- 
thermore, the dissipation of energy within the common envelope is governed by ill-understood 
energy- and angular momentum-transport processes, characterized by a time scale orders of 
magnitude in excess of the local dynamical time scale, and possibly approaching a stellar 
thermal time scale. At present, only the initial immersion of the bin ary components within a 



common envelope can be modeled with a semblance of realism (e.g.. lRicker fc TaamI (120081 )). 



All is not lost, however, in projecting the outcome of common envelope evolution, as the 
entire common envelope process must still satisfy conservation laws for energy and angular 
momentum. Considering first the energy equation, we can write the initial energy of the 
binary, Ei, as 

E, = ^^ + E,, + E2, , (60) 

where Mn is the initial mass of the donor (with envelope intact), A^ is of course the initial 
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binary separation, and En is the initial total energy (gravitational potential plus internal) of 
the donor star, a negative quantity. Subscripts '2i' refer to the corresponding initial quantities 
for the companion star. A similar equation holds for the final, post-common envelope state of 
the binary (subscripts 'f '). As a consequence of common envelope evolution, the envelope of 
the donor star, mass Mcnv = Mu — Mk, is shed from the binary. For reasons outlined in the 
Introduction, we can assume that the companion star survives common envelope virtually 
unchanged in mass or total energy, and so write M2i = M2f = M2 and E2i = E2i = E2. We 
require that the final energy, E^, of the binary be no greater than its initial energy, E-^, which 
yields the inequality 

Af Mif / 2AiA^i \ "^ 

^ < TT^ ( 1 + TTTT^ 1 ' (61) 

where 



A ~ Mii V GMiiMa 



/\Ei = Elf — Eli 

+ U\dm- I + U\dm (62) 

is the difference between final and initial total energies of the donor star, the initial binding 
energy of the ejected envelope. Clearly Mu < Mu, but it is the fact that the initial binding 
energy of the envelope usually greatly exceeds the initial orbital energy of the binary in 
magnitude that is most responsible for dramatic decreases in orbital separation resulting 
from common envelope evolution. 

Considering in turn the conservation of angular momentum, we write the initial angular 
momentum of the binary, Jj , as 



where we explicitly allow for initial orbital eccentricity Cj 7^ 0, but neglect the rotational 
angular momenta of the binary components. The rotational angular momenta are demon- 
strably small, typically contributing less than 2% to the total angular momentum, provided 
the mass ratio of the binary is not too extreme. A parallel expression holds for the final 
angular momentum, Jf. Then requiring, as for energy, that the final angular momentum of 
the binary be no greater than its initial angular momentum, we find 

In general, dissipative processes prior to, or within, common envelope evolution lead us to 
expect that c; ~ ~ Cf. Clearly, save for the possibility that AEi < 0, i.e., that spontaneous 
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ejection by the donor star of its envelope is energetically possible, the energy constraint, 
Eq. fl6T]) . always poses a stronger constraint on post-common envelope orbital separation 
than does Eq. (^^. 

In the years since common envelope evolution became the generally-accepted model for 
the origin of compact bina ries with an evolved component, various approximations have been 
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2000) has been explicitly to integrate over the putative binding energy of the envelope: 
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(65) 



In the limit that this integral extends over the entire mass distribution (Mf — )■ 0), this 
expression gives the total energy of a star of total mass Mj; but it is not rigorously correct 
when parsed over a portion of the star (Mf > 0). The reason is that the total gravitational 
potential energy of a star is a global property of that star. The gravitational potential at 
any point in the star (in the usual convention, where that potential vanishes at infinity) 
depends not only on the mass interior to that point, but also on the mass exterior to it, 
which determines the shape (and integrated depth) of the potential outside that point. 

It may help to clarify the issue to consider an analytic toy model. Consider an isentropic 
n = 1 polytrope. Such a polytrope has the equation of state 



Kp' 



(66) 



where K, the polytropic constant, is a function only of entropy, by assumption. The specific 
internal energy of a gas obeying this equation of state is 



u 



Kp . 



An n = 1 polytrope has the analytic solution (cf. IChandrasekharl Il939l ) 

sin^ 



9 



e 



(67) 



(68) 



where ^ is a dimensionless radial coordinate (^ = r/a), and 6 is the nth (i.e. first) root of the 
ratio of local mass density to central mass density {9 = [p(O/p(0)]^'^")- The surface of the 
polytrope corresponds to the first zero of 9, i.e., ^i = tt. The scale for the radial coordinate 
is, for an n = 1 polytrope, 

K V' 
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and the central density is 



M 



Note (cf. Eq. |69|) that, in the case of an n = 1 polytrope, the radius R = a^i is a function 
only of the polytropic constant K, independent of the mass M. The local mass coordinate 
corresponding to C, is 



M .d^ 

m = ^ — 

71 d^ 

= (4cos4-sm^) 

71 



(71) 



The putative gravitational potential energy per unit mass at radial coordinate ^ is then 

_Gm_GM_ de_ 

r ~ R ^dC ^ ' 

and the internal energy per unit mass 
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where M and R are the total mass and radius of the polytrope, respectively. Integrating 
these energy terms with respect to mass, one finds 
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Thus, at the surface of the polytrope {i = ii = 7i), we find 
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giving for the total energy of the polytrope, 

Etot = n{R) + U{R) = -—— . (78) 

Zri 

Now suppose that all matter outside S, = vr/2 is removed adiabatically from the poly- 
trope. This value of ^ corresponds to mass coordinate m = M/vr and radius coordinate 
r = R/2 in the initial polytrope. Since the specific entropy of the polytrope remains un- 
changed (by assumption), the polytropic constant K, the radial scale factor a, and hence 
the surface radius R all remain unchanged by mass loss. We then have for the final total 
energy of the polytrope 

E,^-^, (79) 

corresponding to a change in total energy 

GM^ f 1 
~2R' V ~^ 



^E = E,-E, = ^^:^(l- — ] . (80) 



But if we integrate over only the envelope that was removed, we find 






2R 



(81) 



that is, the total energy of the initial polytrope! Note that the work done on the polytrope 
in removing this envelope increment by increment from its surface is 

r Gm' , G r,,2 21 

W = - I — — - dm' = — M^ - m^ 
Jm R{m') 2R L J 

i.e., precisely AE, the difference between final and initial total energies of the polytrope. 
Indeed, an adiabatic change of mass is only possible through work done on the polytrope. 

The difference between Eqs. ( IHTj) and ( IHOj) is not particular to this example, but a direct 
result of hydrostatic readjustment in the mass-losing star. As mass is stripped away from 
the surface of the star, a complex readjustment occurs in the interior, as the internal energy 
of the gas does work against the star's self-gravity to re-establish hydrostatic and virial 
equilibrium. Since the surface layers are the coolest part of a star, their loss from the star 
removes little internal energy directly — none at all for a polytropic model, where m — > at 
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the surface. But the loss of matter impacts directly the self-gravitational potential energy 
of the star. Re-establishment of virial equilibrium then converts internal energy from the 
interior into gravitational potential energy, to power expansion to a new state of hydrostatic 
equilibrium. 

The adiabatic mass loss sequences constructed according to the precepts of this paper 
permit us, for the first time, to calculate rigorously the change in total energy between initial 
and final states of the donor star, integrated over their entire interiors. In common envelope 
evolution, this is the change in total energy, AEi, that must be supplied by the dissipation 
of orbital energy and that fixes the largest remnant orbital separation allowed energetically 
(Eq.EH). 



The practical problem of calculating AEi reliably is somewhat more involved. Common 
envelope evolution is of greatest interest, and most prevalent, among binary stars with giant 
branch donors. These stars have very compact cores containing a significant fraction of the 
total stellar mass, and so their cores utterly dominate the evaluation of total energy of these 
stars. In evaluating the change in total energy on loss of the extended giant envelope, one 
is then typically seeking the relatively small difference between large numbers. We estimate 
that our typical relative uncertainty in calculating Ei is of order 2 x 10~^. In comparison, 
the ratio of envelope energy to core energy is (very crudely) of order M'^/R to M^^^^/Rcore] 
with Mcore ~ M/2, and -Rcore ~ 10~^i?, this ratio is 4 x 10~^. Numerical limitations to the 
brute force calculation of the total stellar energy therefore severely limit the utility of that 
approach in just the circumstances when common envelope evolution is of greatest interest. 

Fortunately, the observation above that one can also calculate the change in energy 
of the donor by tracking the energy content of matter removed from its surface provides a 
suitable alternative approach to this problem. Therefore, in addition to calculating AEi by 
direct integration of the total energy of initial and final states (Eq. [^21) . we can calculate the 
change in energy of the donor star by integration along the mass loss sequence: 

\^r + U{R)^ dM , (83) 

where the integrand is evaluated at the surface {r = R, m = M) of the model. This 
expression differs from Eq. (182|) only in that we must account for the finite internal energy 
of matter at the surface of the model; that contribution vanishes at the surface of a classical 
polytrope. Provided we remove small enough mass increments along an adiabatic mass loss 
sequence to track the change in radius smoothly, we achieve a typical relative uncertainty 
in AE[ of order 7 x 10~^. Thus, Eq. fl83l) affords the more accurate measure of the change 
of energy of the donor star early in a mass loss sequence, when \AEi/ Eii\ < 1/40, while 
Eq. fl62l) gives a more accurate measure later in the sequence. In practice, we find the two 
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estimates are in close agreement at the transition, and use a weighted mean of the two to 
determine best estimates of AEi and Ek. 

Given an initial donor of mass Mn and radius Rn, knowledge of AEi{Mif) allows us to 
constrain the range of companion masses, M2, and post-common envelope orbital separations, 
A{, that are energetically allowed. We assume that the initial donor fills its Roche lobe: 

Then the energy constraint (Eq. [6T|) gives the largest post-common envelope orbital sepa- 
ration allowed by conservation of energy, for a post-common envelope (remnant) donor of 

mass Mif! 

, Mif / 2AAE1 Y^ 

To terminate common envelope evolution, we require that the remnant donor fit within its 
Roche lobe, 

and that the companion star also fit within its Roche lobe, 

^^ - t^mJ/M,,) ■ ^''^ 

Fig. [7| illustrates an application of these constraints on survival of a binary with a 3.2 
Mq asymptotic giant branch donor with a 0.80 Mq companion. For that assumed companion 
mass, this figure shows as a function remnant donor mass, Mif the binary separation from 
Eq. (1851) . the binary separation required to accommodate the donor star within its Roche lobe 
from Eq. (1861) . and the corresponding separation required to accommodate the companion 
star within its Roche lobe from Eq. fl87|) . The binary separation must exceed the limits posed 
by the Roche lobe radii for survival of common envelope evolution. Survival of that phase 
does not necessary signal the end of mass transfer, however, only that any further mass 
transfer occur on a thermal time scale or longer. We can ascertain whether the donor likely 
drives continued mass transfer on its thermal time scale upon emergence from the common 
envelope by comparing its nuclear luminosity with its surface luminosity. The difference 
between nuclear and radiated luminosities must be made up work done by expansion or 
contraction within the star; if the nuclear luminosity exceeds the surface luminosity (see 
Fig. Cj), net energy absorption in the stellar interior is expected to drive further thermal time 
scale expansion of the donor. 
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Conclusions and Outlook 



In this paper we have described the construction and interpretation of model sequences 
describing the asymptotic responses of stars to mass loss in a binary system in the adiabatic 
limit, in which mass transfer is so rapid that negligible heat transport occurs within the 
stellar interior during mass loss. We have shown how the radial responses revealed by these 
adiabatic mass loss sequences can be used to establish threshold conditions for dynamical 
time scale mass transfer. That instability may occur promptly upon the donor filling its 
Roche lobe, if the donor has a surface convection zone of any significant depth, or following an 
extended period of thermal time scale mass transfer (a delayed dynamical instability) when a 
radiative donor is stripped down to a nearly isentropic core. Among very luminous stars with 
deep surface convection zones (bright giant branch and asymptotic giant branch stars), the 
adiabatic mass loss sequences show a pronounced initial expansion, as their superadiabatic 
surface layers are stripped away, complicating the interpretation of these models. 

Adiabatic mass loss sequences can also be used to constrain the survival of binaries 
entering common envelope evolution through dynamical time scale mass transfer. We correct 
certain conceptual errors appearing in the literature regarding the calculation of envelope 
binding energy, and describe complementary approaches to its evaluation. By combining 
energy constraints with requirements that both binary components fit within their post- 
common envelope Roche lobes, we can place strict limits on the masses, mass ratios, and 
remnant orbital separations of binaries passing through common envelope evolution. 

In future papers we will survey the Hertzsprung-Russell diagram, mapping out the 
the threshold conditions for dynamical time scale mass transfer, and for the survival of 
common envelope evolution, across the full range of masses and evolutionary stages of 
potential interest, from main sequence stars across the Hertzsprung gap, to giant branch 
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Appendix: Roche Lobe Overflow 

We calculate the mass transfer rate from a lobe-filling star by integrating over a plane 
perpendicular to the line of centers connecting the two binary components, and passing 
through the inner Lagrangian (Li) point. The mass transfer stream pa sses through a sonic 



surface which must coincide closely with this surface of integration (cf. iLubow fc Shulll975 



Eggleton] l2006|). Thus 



Ml = - / pv ■ dS = - / PC dS 

Ju d0 



This preprint was prepared with the AAS I^TJtjX macros v5.2. 
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where p and Cs are the local gas density and sound speed, respectively, and the integral 
is taken over the sonic surface (area element dS), which we assume to coincide with the 
plane of integration defined above. Taylor expansi on of the Roche potential about the inner 
Lagrangian point, Li, point gives, to lowest order ( 1 WebbinkI 1 1 9 771 ) . 



G(Mi + M2) 
A 



--(2a2 + l)(x - xl)' + -(02 - ly 



2«2^ + 



(A2) 



where A is the binary separation; {x,y,z) form a right-handed coordinate system (in units of 
A), centered on the primary (star 1), with the x-axis directed toward the secondary (star 2), 
y-axis in the direction of orbital motion of star 1, and 2;- axis parallel to the orbital angular 
momentum vector; (xl, 0, 0) are the coordinates in this system of Li, and 02 is 



/^ , (l-/i) 



X' 



(A3) 



where /i = Mi /(Mi + M2) = q/{l + q), as before. A simple, accurate approximation to xl is 

(A4) 



XL 



(0.7-0.2gi/3)gi/3 ifg<l, 

1 - (0.7 - 0.2g-i/3j ^-1/3 if g > 1. 



Surfaces of constant (p are, to lowest order, concentric ellpsoids. The intersection of these 
with the plane x = xl are concentric ellipses with axial ratio (02 — l)/a2, and cross-sectional 
area S increasing linearly with potential so that 



2nA^ 



dS 



d(j) G'(Mi + M2) [02(02-1)] 1/2 



(A5) 



We assume for simplicity that the mass outflow from the donor star is laminar, and 
occurs along equipotential surfaces. Let P and p be the pressure and density on a given 
equipotential far from Li, where the material is almost stationary. We assume that they 
follow a power- law adiabat, 

P = Kp^' , (A6) 

where K and Fi = {din P/ din p)^ are constant along streamlines. The n at the sonic su rface 
(x = Xl), the density p^ and sound speed cl on the streamline are (cf. IWebbinkI 119771 ) 



PL 



ri + 1 



i/{ri-i) 



and 



Cl 



or \ 1/2 



Fi + l 



2ri p 



Fi + l p 



1/2 



(A7) 
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Combining Eqs. ( IMj) . flXSJl . (RT]) . and ( IMjl . we can write 

where (ps is the stellar surface potential, and we have used the relation i?L = ^i^A. We have 
combined factors dependent on the mass ratio in the coefficient 

q 1 

r£(l + g) [02(02- l)\^'^ 

a slowly- varying function of q. For the integration over potential 0, we use the approximation 

d0 = ^di?, (All) 

and integrate from _Rl to the stellar radius R. 

Eq. (]A9P is at best a rough approximation. In addition to approximation of the equation 
of state by a simple power law, neither of the assumption of laminar flow nor that of flow 
along equipotential surfaces is valid. A substantial fraction of the mass outflow well s up from 



within the Roche lobe of the donor, rather than flowing along its surface (see, e.g., [Eggleton 



( 120061 )). However, we find that the adiabatic index Fi often varies dramatically with depth 
through the outer envelope of the donor star, and when Fi < 1 (as often happens in partial 
ionization zones) no solution for an irrotational outfiow exists. The outfiow must then be 
turbulent. The assumptions of laminar fiow and equipotential fiow err in opposite senses. 
A significant fraction of the work done by pressure gradients is spent in driving turbulence 
rather than bulk fiow, reducing the mass transfer rate, while the upwelling of gas from the 
interior increases it. We consider that our estimate for Mi above is accurate within a factor 
of two or so — for constant F i > 1, it yields the same functional dependence on Mi, i?L, 



and 0s — 0L as does lEggletonl 's (120061) prescription — but a rigorous treatment will require 



full three-dimensional models of the the outfiow. 
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Fig. 1. — Specific entropy profiles at five epochs in a 2.0 Mq donor undergoing conservative 
Roche lobe overfiow (RLOF) to a 0.65 M© companion. Curve a marks the beginning of 
RLOF. The mass-transfer reaches rate quickly rises from 5.5 x lO^^MQyr^^, roughly a 
thermal time scale rate (curve 6), to 5.1 x 10^^ M^gyr"-*^ (curve e). 
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Fig. 2. — Density (g cm ^), luminosity, temperature (K) and radius profiles for the 2.0 Mq 
donor shown in Fig. [T] at epochs c (sohd lines) and e (dotted lines). 
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Fig. 3. — Adiabatic response curves (top) and the entropy profiles (bottom) for low- mass 
ZAMS models, of initial mass M^ = 0.30, 0.40, . . ., 0.90, 1.00 Mq. Solid circles mark the 
boundaries between the radiative core and the convective envelope. Initial masses increase 
from top to bottom in the radius response curves (upper panel), and from bottom to top in 
the entropy profiles (bottom panel). 
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Fig. 4. — Radial response curves for a 3.2 Mq terminal main sequence donor star. The solid 
line traces the adiabatic mass loss sequence. Broken lines mark the Roche lobe radius as 
a function of mass for three initial mass ratios: go = ^ii/^2i = 3.0 (dotted line), 3.757 
(dash-dotted line, and 4.5 (dashed line), all assuming conservative mass transfer. In the 
limit that the donor star truly responds to mass loss adiabatically, the donor star becomes 
unstable to dynamical time scale mass transfer when the adiabatic radius exceeds the Roche 
lobe radius. Thus, the go = 3.0 binary is stable against dynamical mass transfer throughout 
mass loss, while the go = 4.5 binary reaches dynamical instability when the donor star has 
been stripped to Mi/Mn = 0.853. This is a delayed dynamical instability, because it erupts 
only after an extended period of thermal time scale mass tranfer. The threshold initial mass 
ratio for delayed dynamical instability occurs where the Roche lobe radius and adiabatic 
radius reach tangency, roughly at Mi/Mn = 0.653, go = 3.757. 
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Fig. 5. — Mass loss rates of the 3.2Mq terminal main sequence models of Fig. |4]from fully 
time-depen dent evolut i onary calculations of conservative Roche lobe overflow (based on the 
method by iHan et al.l (l2000f )). The onset of dynamical time scale mass transfer from these 
calculations occurs at Mi/Mu < 0.905 for go = 4.5 and Mi/Mu < 0.814 for go = 3.757. 
These calculations stop short of fully djTiamical time scale mass transfer, but suggest that 
the time-dependent models may be slightly more prone to instability. 
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Fig. 6. — Adiabatic mass loss sequences for true (thick solid lines) and pseudo-models (thick 
dashed lines) of a 1 Mq star near the tip of the giant branch. The pseudo-models have 
isentropic envelopes, enabling us to quantify Aexp, the prompt initial expansion resulting 
from removal of superadiabatic surface layers. Also shown are the Roche lobe radii within 
true and pseudo-models (thin solid and dashed lines, respectively) at which the mass loss 
rate would reach Mkh- At higher mass loss rates, thermal relaxation freezes out. Critical 
initial mass ratio at which Mi just reaches Mkh for the pseudo- model is shown by the heavy 
dotted. (The corresponding limit for the true model is omitted for clarity). We believe 
that these two measures of the critical mass ratio probably bracket the true threshold for 
dynamical instability. 
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Fig. 7. — Limiting remnants of common envelope evolution. In this example, the donor is 
a 3.2 Mq asymptotic-branch donor star, and its companion a 0.8 M© main sequence star. 
The solid line marks the maximum post-common-envelope separations allowed by energy 
conservation (Eq. |85|) . the dashed line the minimum separation needed to accommodate the 
donor within its Roche lobe (Eq. [861 ^ind the dotted line the minimum separation needed 
to accommodate the companion star. No systems with M{ > 0.690 M© survive, because 
the donor star radius is still too large to fit within its Roche lobe, while systems with 
Mf < 0.663 Mq are too compact to accommodate the companion star. Note also that 
remnant donors with with masses Mk < 0.75 Mq have nuclear luminosities exceeding their 
radiated luminosities; they must expand beyond the dashed line to reach thermal equilibrium, 
and so may drive a further phase of thermal time scale mass transfer. 



